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We assemble a sample of 3258 low-redshift galaxies from the Sloan Digital Sky Survey Data 
Release 6 (SDSS DR6) with complementary photometric observations by GALEX, 2MASS 
and IRAS at far-ultraviolet and infrared wavelengths. We use a recent, simple but physically 
C$ \ motivated model to interpret the observed spectral energy distributions of the galaxies in this 

sample in terms of statistical constraints on physical parameters describing the star formation 
history and dust content. The focus on a subsample of 1658 galaxies with highest signal-to- 
noise ratio (S/N) observations enables us to investigate most clearly several strong correlations 
between various derived physical properties of galaxies. In particular, we find that the typical 
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f- ») ' dust mass M<j of a galaxy forming stars at a rate tp can be estimated remarkably well using the 

formula M d = (1.28 ± 0.02) x 10 7 (ip/M® yr -i)i-U±o.oi Mq over at least three orders of 
{Sj ■ magnitude in both quantities. We also find that the dust-to-stellar mass ratio, the ratio of dust 

mass to star formation rate and the fraction of dust luminosity contributed by the diffuse in- 
terstellar medium (ISM) all correlate strongly with specific star formation rate. A comparison 
with recent models of chemical and dust evolution of galaxies suggests that these correlations 
could arise, at least in part, from an evolutionary sequence. As galaxies form stars, their ISM 
becomes enriched in dust, while the drop in gas supply makes the specific star formation rate 
decrease. Interestingly, as a result, a young, actively star-forming galaxy with low dust-to-gas 
k> , ratio may still be highly dusty (in the sense of a high dust-to-stellar mass ratio) because it 

'j_J ■ contains large amounts of interstellar gas. This may be important for the interpretation of the 

infrared emission from young, gas-rich star-forming galaxies at high redshift. The results pre- 
sented in this paper should be especially useful to improve the treatment of the ISM properties 
of galaxies in semi-analytic models of galaxy formation. Our study also provides a useful lo- 
cal reference for future statistical studies of the star formation and dust properties of galaxies 
at high redshifts. 

Key words: dust, extinction - galaxies: ISM - galaxies: stellar content - galaxies: statistics 
- galaxies: evolution. 



1 INTRODUCTION 

Interstellar dust and star formation activity are strongly linked in 
galaxies. Dust grains condense in the cold envelopes of evolved 
stars and in supernova ejecta. In return, they favour the formation 
of molecular hydrogen, shield the newly-formed molecules from 
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ultraviolet radiation and participate in the formation and cooling of 
molecular clouds, which collapse to form new stars. 

Studies of the ultraviolet, optical and near-infrared emission 
from large samples of local galaxies have shed some light on the re- 
lation between star formation activity and dust content. Such stud- 
ies show that, in general, galaxies with the highest star formation 
rates also suffer the highest ultraviolet and optical attenuation (e.g., 
Kauffmann et al. 2003b; Brinchmann et al. 2004). However, obser- 
vations at ultraviolet, optical and near-infrared wavelengths set only 
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limited constraints on the dust content of galaxies. New insight into 
the dust properties requires observations at mid- and far-infrared 
wavelengths, where the dust re-radiates the energy absorbed from 
starlight. The first studies of this kind combined ultraviolet and op- 
tical observations of large samples of nearby galaxies with empiri- 
cal estimates of the total infrared luminosity computed from obser- 
vations with the Infrared Astronomical Satellite (IRAS) at 12, 25, 
60 and 100 /ttm (e.g. Wang & Heckman 1996; Hopkins et al. 2001; 
Sullivan et al. 2001; Kong et al. 2004). Other recent studies have 
taken this (largely empirical) approach a step further by develop- 
ing consistent models of galactic spectral energy distributions at 
ultraviolet, optical and infrared wavelengths (e.g. Silva et al. 1998; 
Dopita et al. 2005). These sophisticated models combine the spec- 
tral evolution of stellar populations with detailed calculations of 
the transfer of starlight through the interstellar medium in galaxies. 
Such models, however, are not optimised to interpret the observa- 
tions of large samples of galaxies. 

Recently, da Cunha et al. (2008) proposed a simple, physi- 
cally motivated model to interpret consistently the ultraviolet, opti- 
cal and infrared observations of large samples of galaxies in terms 
of statistical constraints on physical parameters pertaining to the 
stars, gas and dust. By design, this model is less sophisticated 
than other existing radiative transfer models. In return, it allows 
one to interpret efficiently the spectral energy distributions of large 
samples of galaxies in terms of statistical constraints on a mini- 
mal set of adjustable physical parameters, da Cunha et al. (2008) 
used this model to derive median-likelihood estimates of the stel- 
lar mass, star formation rate, dust attenuation and dust mass from 
the observed ultraviolet, optical and infrared spectral energy distri- 
butions of 66 galaxies in the Spitzer Infrared Nearby Galaxy Sam- 
ple (SINGS, Kennicutt et al. 2003). For these galaxies, the dust- 
to-stellar mass ratio appears to correlate strongly with specific star 
formation rate (i.e. the star formation rate divided by the total stellar 
mass), confirming the expectation that dust mass and star formation 
rate are tightly related in galaxies. 

In the present paper, we investigate further the relation be- 
tween star formation activity and dust content by studying the prop- 
erties of a much larger sample of 3258 star-forming galaxies, for 
which photometric observations are available at ultraviolet, opti- 
cal and infrared wavelengths. We use the model of da Cunha et al. 
(2008) to derive median-likelihood estimates of several physical 
parameters of each galaxy in this sample, such as the total stel- 
lar mass, dust mass star formation rate and fractional contribu- 
tion of different dust components to the total infrared luminos- 
ity. We then focus on a sub-sample of 1658 galaxies with high- 
est signal-to-noise ratio (S/N) photometry from the far-ultraviolet 
to the far-infrared. The large size of this sample allows us to 
study in unprecedented detail the relations between specific star 
formation rate, dust mass and dust-to-gas ratio already identified 
in the SINGS sample by da Cunha et al. (2008). We find that, 
in particular, the typical dust mass Md of a galaxy forming stars 
at a rate tp can be remarkably well estimated using the formula 
M d = (1.28 ± 0.02) x 10 7 O/M yr -i)i n±o oi Mq _ 

We interpret our results in the framework of recent models of 
chemical and dust evolution of galaxies by Calura et al. (2008). 
These models provide a means of predicting the time evolution of 
the stellar and dust content of galaxies with different star formation 
histories. We conclude from this that the relations between specific 
star formation rate and dust content exhibited by the galaxies in 
our sample could be regarded, at least in part, as an evolutionary 
sequence. 

The paper is organised as follows. We first describe the 



matched GALEX-SDSS-2MASS-IRAS sample in Section 2. In Sec- 
tion 3, we outline the method used to derived statistical estimates of 
the total stellar mass, dust mass, star formation rate and fractional 
contribution of different dust components to the total infrared lu- 
minosity for each galaxy in this sample. We present the relations 
between these various derived quantities in Section 4. In Section 5, 
we discuss some potential biases of our approach. We also compare 
our results with the predictions of different models to illustrate the 
implications of our study for the evolution of the dust content of 
star-forming galaxies. Section 6 summarises our conclusions. 



2 THE GALAXY SAMPLE 

In this work we analyse a galaxy sample obtained by cross- 
correlating the Sloan Digital Sky Survey Data Release 6 (SDSS 
DR6) main spectroscopic sample with photometric catalogues at 
ultraviolet (from the Galaxy Evolution Explorer, GALEX), near- 
infrared (from the Two Micron All Sky Survey, 2MASS) and far- 
infrared (from the Infrared Astronomical Satellite, IRAS) wave- 
lengths. 



2.1 Optical photometry 

The SDSS DR6 spectroscopic sample (Adelman-McCarthy et al. 
2008) contains 792,680 galaxies with known redshifts to a Pet- 
rosian r-band magnitude limit r < 17.77. Optical photometry in 
the ugriz bands (at 3557, 4825, 6261, 7672 and 9097 A, respec- 
tively) is available for these galaxies. We use the SDSS 'model 
magnitudes' in these bands, which reflect the integrated light from 
the whole galaxy and are the best suited to comparisons with 
total photometry from GALEX, 2MASS and IRAS at ultraviolet 
and infrared wavelengths. We stress that such consistent multi- 
wavelength photometry is a key requirement in the interpretation 
of the spectral energy distributions of the galaxies in Section 3. 

Many fundamental parameters of the galaxies of the SDSS 
DR6 spectroscopic sample have been derived from the spectro- 
scopic data (e.g. Brinchmann et al. 2004), and are publicly avail- 
able in an online database 1 . We use this information to restrict our 
sample to galaxies classified as 'star-forming' according to their 
emission lines. Active Galactic Nuclei (AGNs) in the sample are 
eliminated when at least one of two criteria are met: (i) the presence 
of broad emission lines from high-velocity gas around the AGN 
(e.g. Brinchmann et al. 2004); (ii) the presence of narrow emission 
lines with intensity ratios characteristic of an AGN in the Baldwin, 
Phillips & Terlevich diagram (BPT; Baldwin et al. 1981), accord- 
ing to the criterion of Kauffmann et al. (2003a) 2 . If the emission 
lines are highly attenuated, the BPT diagnostic may fail to identify 
an obscured AGN (in Section 5.3, we discuss the possible contam- 
ination of our sample by optically thick AGN). Observations in the 
X-rays or in the mid-infrared can be used to detect AGN in this 
case. 



1 http://pc-66.astro.up.pt/~jarle/GASS/ 

2 The BPT diagram allows one to distinguish star-forming galaxies from 
AGN on the basis of the degree of ionisation of the gas by comparing the 
narrow-line ratios [O III]A 5007 A/H/3 against [N II]A 6584 A/H«. 
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2.2 Ultraviolet photometry 

To supplement SDSS DR6 spectroscopic data with ultraviolet mea- 
surements, we cross-correlate the sample with the latest data re- 
lease from GALEX (Martin et al. 2005; Morrissey et al. 2005). 
GALEX is an all-sky survey providing images of the galaxies in 
two photometric bands: the far-ultraviolet (FUV, at 1520 A) and 
the near-ultraviolet (NUV, at 2310 A). To match the SDSS and 
GALEX catalogues, we use a search radius of 4 arcsec around the 
SDSS position (data kindly provided in advance of publication by 
David Schiminovich). We require a detection in at least one of the 
two FUV and NUV bands. We retain galaxies which have a single 
GALEX detection within the search radius, i.e., we eliminate ob- 
jects for which two or more GALEX detections may exist for the 
same SDSS source. We adopt GALEX FUV and NUV 'automags', 
which use a Kron aperture defined by the profile of each galaxy in 
each band, and are designed to recover the total flux. 

At ultraviolet wavelengths, Galactic foreground extinction is 
particularly important. Therefore, we apply Galactic reddening cor- 
rections to the GALEX magnitudes (Seibert et al. 2005): Apuv = 
8.29 E(B - V) and A MV = 8.18 E(B - V), where the colour 
excess E(B — V) is determined from the dust reddening maps 
of Schlegel et al. (1998). Although observations at longer wave- 
lengths are less affected by Galactic extinction, we also correct the 
observed optical magnitudes of Section 2. 1 and near-infrared mag- 
nitudes of Section 2.3 using the same method. 
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Figure 1. Total infrared luminosity L ^ ot from equation ( 1 ) plotted against 
redshift for galaxies in the matched GALEC-SDSS-2MASS-/SAS described 
in Section 2. Galaxies belonging to the IRAS PSCz catalogue (blue sym- 
bols) are distinguished from those belonging to the IRAS FSC catalogue 
(red symbols). 



2.3 Near-infrared photometry 

We further supplement our sample with near-infrared photometry 
in the JHK S bands (at 1.25, 1.65 and 2.17 /im, respectively) from 
the 2MASS (Skrutskie et al. 2006) All Sky Extended Source Cata- 
log (XSC). We cross-identify the SDSS DR6 spectroscopic sample 
and the 2MASS XSC within a search radius of 5 arcsecs around 
the SDSS coordinates. Only sources with no artefacts and which 
are not in close proximity to a large nearby galaxy are retained. 
Concerning the photometry, we follow the recommendations in the 
User's Guide to the 2MASS All-Sky Data Release 3 to include most 
of the integrated flux from the galaxies while providing accurate 
colours: we adopt K20 fiducial isophotal elliptical aperture mag- 
nitudes for galaxies with K 3 < 13, and fixed 7 arcsec circular 
aperture magnitudes for fainter galaxies. 

2.4 Mid- and far-infrared photometry 

We now turn to the most original feature of our sample: the in- 
clusion of mid- and far-infrared photometry obtained with IRAS 
(Beichman et al. 1988). Observations with IRAS are rather limited 
in sensitivity, but they have the advantage of providing an all-sky 
survey at 12, 25, 60 and 100 ^m. 4 We cross-correlate our sam- 
ple with both the IRAS Point Source Catalogue (PSCz; Saunders et 
al. 2000) and the IRAS Faint Source Catalogue v2.0 (FSC; Moshir 
1989). The PSCz catalogue is a complete and uniform galaxy cata- 
logue assembled from the IRAS Point Source Catalogue (Beichman 
et al. 1988) and supplemented by various redshift surveys. It con- 
tains 15,411 galaxies with measured redshifts to a depth of 0.6 Jy 
at 60 /im. The Faint Source Catalogue contains 173,044 sources to 

3 http://www.ipac.caltech.edu/2mass/releases/allsky/doc/ 

4 In practice, most measurements at 12 and 25 /im are upper limits on the 
flux density, and we do not include them in our spectral fits. Only about 3 
per cent of the galaxies in the sample are detected in all 4 IRAS bands. 



a depth of roughly 0.25 Jy at 60 (im. Low Galactic latitude regions 
(|6| < 20°) are excluded from this catalogue because of the con- 
tamination by foreground Galactic sources at this detection limit. 
For both catalogues, we rely on the quoted data quality flags to 
select sources having a good flux quality at 60 /j,m and at least a 
moderate flux quality at 100 /im; detections at 12 and 25 /im are 
not strictly required. To increase the reliability of the sample, we 
exclude observations which are confusion-limited or contaminated 
by cirrus emission. 

We cross-correlate the SDSS sample with the IRAS PSCz 
and IRAS FSC catalogues using the same criteria as Pasquali 
et al. (2005), i.e., we discard IRAS sources with more than one 
SDSS-matched object within a search radius of 30 arcsec. For the 
PSCz catalogue (with redshift information), this criterion appears 
to cause only about 1 per cent of the infrared sources to be assigned 
the wrong SDSS counterpart. The contamination may be higher in 
the case of the lower-quality FSC catalogue. Using simulated data, 
Pasquali et al. (2005) estimate that the percentage of wrong cross- 
identifications is at most 1.5 per cent for a matching radius of 30 
arcsec. 

2.5 Final sample 

Our final GALEX-SVSS-2MASS-IRAS sample is composed of 
3258 galaxies from the PSCz and FSC catalogues at redshifts 
z ^ 0.22 (see Fig. 2a). To summarise, the galaxies in this sam- 
ple have: 

• ultraviolet fluxes in at least one of the FUV and NUV GALEX 
bands; 

• optical ugriz fluxes from SDSS; 

• near-infrared JHK S fluxes from 2MASS; 

• 60- and 100- flux densities from IRAS (some galaxies also 
have 12- and 25-/im measurements). 
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We note that, for consistency, all fluxes used in our study are de- 
signed to be total fluxes, which do not require aperture correc- 
tions. For this reason, we also choose not to include SDSS spectro- 
scopic information available for these galaxies, which is restricted 
to the inner 3-arcsec diameter aperture sampled by the SDSS spec- 
troscopic fibre aperture. At the low redshifts of our galaxies, this 
limited spatial sampling could severely bias estimates of the star 
formation rate and dust content (Kewley et al. 2005). In addition, 
spectroscopic quantities pertaining to a restricted central area can- 
not be compared directly with multi-wavelength photometric quan- 
tities describing the whole galaxy. 

It is instructive to examine the typical infrared luminosity L\ ot 
of the galaxies in our sample. For this purpose, we compute L^ ot 
from the 60- and 100-/xm IRAS flux densities and Lj 00 using 
the empirical formula (Helou et al. 1988) 

Lj ot = F c Lfir , (1) 
where 

L FIR = 1.26 x 10~ 14 (2.58 F„ 60 + L„ 100 ) 4ndt . (2) 

Here, di. is the luminosity distance in m, Ft and Fl 00 are in Jy, 
and Lj ot and Lfir are in W. We compute the correction factor 
F c to obtain L^ ot from Lfir in equation (1) following the em- 
pirical prescription of Helou et al. (1988), which depends on the 
F„° /F^ 00 ratio, and assuming a dust emissivity index ft = 2. The 
resulting median correction for our sample is F c = 1.35. 

In Fig. 1, we plot the total dust infrared luminosity L\ ot com- 
puted in this way as a function of redshift. A large fraction of the 
galaxies in our sample (about 22 per cent) have L\ ot > 10 11 Lq. 
Such galaxies are usually referred to as 'luminous infrared galax- 
ies' (LIRGs). About 1 per cent of the galaxies of our sample have 
Lj ot > 10 12 Lq. These galaxies with extremely high infrared lu- 
minosities are usually referred to as 'ultra-luminous infrared galax- 
ies' (ULIRGs). This type of galaxies has been the object of exten- 
sive studies (e.g. Soifer et al. 1987; Veilleux et al. 1995, 1999; 
Rigopoulou et al. 1999; Cao et al. 2006; Armus et al. 2007). In 
the local universe, most ULIRGs are observed to be the results of 
interactions and mergers. Visual inspection of the SDSS optical im- 
ages confirms that ULIRGS in our sample also have disturbed mor- 
phologies. 

In Fig. 2, we compare the properties of our sample with the 
overall properties of the SDSS star-forming galaxy sample. Fig. 2a 
shows that galaxies in our sample tend to lie at lower redshifts than 
the bulk of SDSS star-forming galaxies. This results from the re- 
quired detection by low-sensitivity IRAS observations. Fig. 2b fur- 
ther shows that galaxies in our sample tend to have typically red- 
der g — r colours than SDSS star-forming galaxies. This is most 
probably a consequence of the higher dust content of our galaxies, 
which makes them detectable by IRAS. In fact, Obric et al. (2006, 
see their figure 19) point out that SDSS galaxies detected by IRAS 
have systematically higher dust attenuation than the average SDSS 
star-forming galaxy. The difference in the distributions of r-band 
absolute magnitude M r in Fig. 2c is more subtle to interpret. This 
is illustrated by Fig. 2d, where we plot M r in different redshift bins 
for both samples. At redshifts z < 0.12, galaxies in our sample are 
typically brighter than SDSS star-forming galaxies. This is because 
SDSS galaxies with intrinsically faint r-band magnitudes are too 
faint in the infrared to be detected by IRAS. At redshifts z > 0.12, 
the SDSS and IRAS detection limits both correspond to similarly 
bright galaxies. 




Figure 2. Properties of the SDSS DR6 star-forming galaxies (in grey) 
and the matched GALEX-SDSS-2MASS-IRAS subsample considered in this 
work (in black). The histograms show the normalised distributions of the 
following parameters: (a) Redshift, z; (b) Galaxy g — r colour, fc-corrected 
to redshift z = 0.1; (c) Absolute r-band model magnitude, fc-corrected to 
redshift z = 0.1, Mj?- 1 . Panel (d) shows the mean r-band absolute magni- 
tude, M r , in different redshift bins for both samples (the error bars represent 
the standard deviation ±<7 in each bin). 



3 STATISTICAL CONSTRAINTS ON PHYSICAL 
PARAMETERS 

In this Section, we use the simple model of da Cunha et al. (2008) 
to extract star formation histories and dust contents from the ul- 
traviolet, optical and infrared observations of the galaxy sample 
described in Section 2. We first briefly summarise the model. 
Then, we describe the statistical approach used to derive median- 
likelihood estimates of physical parameters from the available data. 

3.1 Description of the model 

The simple, physically motivated model of da Cunha et al. (2008) 
allows us to interpret the mid- and far-infrared spectral energy dis- 
tributions of galaxies consistently with the emission at ultraviolet, 
optical and near-infrared wavelengths. We briefly recall the main 
features of this model. 

We compute the emission by stars in galaxies using the lat- 
est version of the Bruzual & Chariot (2003) population synthe- 
sis code (Chariot & Bruzual, in preparation). This code predicts 
the spectral evolution of stellar populations in galaxies from far- 
ultraviolet to far-infrared wavelengths and at ages between 1 x 10 5 
and 2 x 10 10 yr, for different metallicities, initial mass functions 
(IMFs) and star formation histories. In this work, we adopt the 
Chabrier (2003) Galactic-disc IMF. 

The emission from stars is attenuated using the simple two- 
component dust model of Chariot & Fall (2000). This accounts for 
the fact that stars are born in dense molecular clouds with typical 
lifetimes of 10 7 yr; at later ages, stars migrate to the ambient (dif- 
fuse) ISM. Thus, the light produced by stars younger than 10 7 yr 
is attenuated by dust in the birth clouds and in the ambient ISM, 
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while the light produced by older stars is attenuated only by dust 
in the ambient ISM. The model of Chariot & Fall (2000) uses an 
'effective absorption' curve for each component, t\ oc X~ n . The 
slope n reflects both the optical properties and the spatial distribu- 
tion of the dust. Following Chariot & Fall (2000), we adopt for the 
ambient ISM 



- ISM 



^tV(A/5500 A)" 



(3) 



where fy is the total effective F-band absorption optical depth 
of the dust seen by young stars inside birth clouds, and /i — 
Tv SM /{tv C + Ty SM ) is the fraction of this contributed by dust 
in the ambient ISM. For the stellar birth clouds, we adopt: 

f A BC = (1 - n)fv (A/5500 A)" 1 ' 3 . (4) 

We use this prescription to compute the total energy absorbed by 
dust in the birth clouds and in the ambient ISM; this energy is re- 
radiated by dust at infrared wavelengths. By analogy with da Cunha 
et al. (2008), we define the total dust luminosity re-radiated by dust 
in the birth clouds and in the ambient ISM as L d ° and Lj SM , 
respectively. The total luminosity emitted by dust in the galaxy is 
then 



T BC . T ISM 
J-'d T J->d 



(5) 

We distribute L d c and Lj SM in wavelength over the range 
from 3 to 1000 /im using 4 main components (see da Cunha et al. 
2008 for detail): 

• the emission from polycyclic aromatic hydrocarbons (PAHs; 
i.e. mid-infrared emission features), 

• the mid-infrared continuum emission from hot dust with tem- 
peratures in the range 130-250 K, 

• the emission from warm dust in thermal equilibrium with ad- 
justable temperature in the range 30-60 K, 

• the emission from cold dust in thermal equilibrium with ad- 
justable temperature in the range 15-25 K. 

In stellar birth clouds, the relative contributions to L d ° by PAHs, 
the hot mid-infrared continuum and warm dust are kept as ad- 
justable parameters. These clouds are assumed not to contain any 
cold dust. In the ambient ISM, the contribution to L d SM by cold 
dust is kept as an adjustable parameter. The relative proportions of 
the other 3 components are fixed to the values reproducing the mid- 
infrared cirrus emission of the Milky Way. da Cunha et al. (2008) 
find that the above minimum number of components is required to 
account for the infrared spectral energy distributions of galaxies in 
a wide range of star formation histories. 



3.2 Median-likelihood estimates of physical parameters 

The model summarised in Section 3.1 allows us to derive statistical 
estimates of galaxy physical parameters, such as the stellar mass, 
star formation rate and dust mass, from simultaneous fits of ultra- 
violet, optical and infrared observations. To achieve this, we adopt 
a Bayesian approach similar to that used by da Cunha et al. (2008) 
to interpret the spectral energy distributions of SINGS galaxies. 



3.2.1 Model library 

We build large libraries of stochastic models at different redshifts 
z = 0.00, 0.05, 0.10, 0.15 and 0.20. 

At each redshift, we generate a random library of stellar popu- 
lation models for wide ranges of star formation histories, metallici- 
ties and dust contents. Each star formation history is parameterised 



in terms of an underlying continuous model with exponentially de- 
clining star formation rate, on top of which are superimposed ran- 
dom bursts (see also Kauffmann et al. 2003b). The models are dis- 
tributed uniformly in metallicity between 0.2 and 2 times solar. The 
attenuation by dust is randomly sampled by drawing the total effec- 
tive F-band absorption optical depth, fv, between and 6, and the 
fraction of this contributed by dust in the ambient ISM, /j,, between 
and 1 (see da Cunha et al. 2008 for more details on the prior dis- 
tributions of these parameters). For each model, we compute the 
fraction of the total energy absorbed by dust in the diffuse ISM, / M . 
We also compute the specific star formation rate averaged over the 
last t 8 = 10 s yr, 



Mt) = 



t 8 M„{t) 



(6) 



where M* (t) is the stellar mass at time t. 

In parallel, at each redshift, we generate a random library of 
infrared spectra for wide ranges of dust temperatures and fractional 
contributions by different dust components to the total infrared lu- 
minosity. The fraction of total dust luminosity contributed by the 
diffuse ISM, / M = L] SM /L d ot , and the fractional contribution 
by warm dust in thermal equilibrium to the total dust luminos- 
ity of the birth clouds, £w C > are uniformly distributed between 
and 1 ; the fractions of luminosity contributed by the other compo- 
nents in the birth clouds, £mir ( n °t mid-infrared continuum) and 
£f?yh (PAHs), are randomly drawn from prior distributions such 
that £w C + £mir + £pah = 1- We distribute uniformly the temper- 
ature of warm dust in thermal equilibrium in the birth clouds, T^°, 
between 30 and 60 K, and the temperature of cold dust in thermal 
equilibrium in the diffuse ISM, Tc SM , between 15 and 25 K. Fi- 
nally, we distribute uniformly the contribution by cold dust to the 
total infrared luminosity of the ISM, £<i SM , between 0.5 and 1. 

We compute the dust mass associated to each model in the 
library of infrared spectra as 



M d 



1.1 (M BC 



M^ M + M<f M ) , 



(V) 



where , M™ and Mq SM are the masses of the dust com- 
ponents in thermal equilibrium (warm dust in the birth clouds and 
the ambient ISM and cold dust in the ambient ISM). The multiply- 
ing factor 1 . 1 accounts for the small contribution by stochastically 
heated dust (see da Cunha et al. 2008 for detail). The mass M d (T d ) 
of dust in thermal equilibrium at the temperature T d is estimated 
from the corresponding far-infrared luminosity, L^ d , using the re- 
lation (Hildebrand 1983) 



L A Td = An M d (T d ) k x B x (T d ) 



(8) 



where k\ is the dust mass absorption coefficient and B\(T d ) the 
Planck function of temperature T d . Following da Cunha et al. 
(2008), we adopt k,\ oc A - * 3 , with /3 = 1.5 for warm dust and 
P — 2.0 for cold dust, normalised to Kgso = 0.77 g _1 cm 2 at 
850 /jm (Dunne et al. 2000). 

We combine the stochastic libraries of attenuated stellar spec- 
tra and dust emission spectra at each redshift by associating models 
corresponding to similar values of / M (to within some uncertainty 
interval <5/ M = 0.15) in the two libraries, which we scale to the 
same total dust luminosity L d ot . For each combined spectrum, we 
compute the synthetic photometry in the GALEX FUV and NUV, 
SDSS ugriz, 2MASS JHK S and IRAS 12-, 25-, 60- and lOO-^m 
bands. 
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3.2.2 Corrections applied to the observed fluxes 

In the redshift range of the sample described in Section 2, promi- 
nent optical nebular emission lines can significantly affect the ob- 
served galaxy fluxes in the SDSS gri bands. The model spectra we 
use to interpret these data do not include nebular emission lines. 
Therefore, to interpret the optical fluxes of observed galaxies with 
these models, we first need to correct the observed SDSS gri mag- 
nitudes for potential contamination by nebular emission lines (e.g., 
Kauffmann et al. 2003b). We use the corrections inferred by Jarle 
Brinchmann (private communication) from fits of the stellar con- 
tinuum emission of each SDSS optical spectrum with Bruzual & 
Chariot (2003) models (we assume for simplicity that the correc- 
tion derived in this way within the aperture sampled by the SDSS 
fibre applies to the galaxy as a whole). 

We also compute ^-corrections to the ultraviolet, optical and 
near-infrared magnitudes of each galaxy in the sample. To min- 
imise these, we fc-correct the magnitudes from the galaxy redshift 
to the closest redshift of the model grid described in Section 3.2.1, 
i.e., z = 0.00, 0.05, 0.10, 0.15 or 0.20 (we use the version v3 of 
the KCORRECT code of Blanton et al. 2003). This procedure cannot 
be extended to the IRAS 12-, 25-, 60- and 100-^im flux densities, to 
which we do not apply any ^-correction. We do not expect this to 
have any noticeable influence on our results, given the large effec- 
tive width of the IRAS filter response functions and the relatively 
large observational uncertainties in these bands. 

To account for the uncertainties linked to the fc-correction and 
emission-line correction, we add the following errors to the quoted 
flux uncertainties: 2 per cent for GALEX, 2MASS, and SDSS z 
bands, and 1.5 per cent for the SDSS gri bands. For the less ac- 
curate SDSS u-band photometry, we take an overall observational 
uncertainty of 10 per cent. 



3.2.3 Spectral fits 

We perform spectral fits by comparing the observed spectral en- 
ergy distribution of a galaxy to every model in the library at the 
corresponding redshift. Specifically, for each observed galaxy, we 
compute the \ 2 goodness of fit of each model. A model is charac- 
terised by a set of randomly drawn physical parameters. We build 
the likelihood distribution of any given physical parameter for the 
observed galaxy by weighting the value of that parameter in each 
model by the probability cxp(— x 2 /2). We take our final estimate 
of the parameter to be the median of the likelihood distribution, and 
the associated confidence interval to be the 16th-84th percentile 
range. 

We use this approach to derive the likelihood distributions of 
several physical parameters of the galaxies in our sample, based on 
fits of the GALEX FUV and NUV, SDSS ugriz, 2MASS JHK„ 
and IRAS 12-, 25-, 60- and 100-/im fluxes. We focus particularly 
on: the star formation rate averaged over the last 10 8 yr, ip; the 
stellar mass, M„; the specific star formation rate, ips = ip/M t ; 
the dust mass, Md\ the total luminosity of the dust, L^ ot , and the 
fraction of this contributed by the diffuse ISM, / M . 

We first check how well the model can reproduce the observed 
spectral energy distributions of the galaxies in our sample. The his- 
tograms in Fig. 3 show, for each photometric band, the distribu- 
tion of the difference between the observed luminosity L° bB and 
the best-fit model luminosity £™ od , in units of observational er- 
ror a. Overall, the model provides remarkably consistent fits to the 
observed ultraviolet, optical and infrared luminosities of the galax- 
ies. Fig. 3 shows small systematic offsets in the g, r and z bands, 




Figure 3. Distribution of the difference between observed luminosity L° 
and best-fit model luminosity L™ od , in units of the observational error a, 
for the galaxies in the matched GAL£X-SDSS-2MASS-ffiAS sample de- 
scribed in Section 2. Each panel refers to a different photometric band, as 
indicated. For reference, the black dashed curve shows a Gaussian distribu- 
tion with unit standard deviation. The best-fit model for each galaxy was se- 
lected by fitting as many fluxes as available in the following bands: GALEX 
(FUV and NUV), SDSS (ugriz), 2MASS (JHK S ), and IRAS (12, 25, 60 
and 100 /an). 



corresponding to an overestimate of the p-band flux and underesti- 
mates of the r- and 2-band fluxes of the order of 0.01 mag. These 
offsets may originate from a deficiency in stellar population syn- 
thesis models, further worsened by the potential contamination of 
the g and r bands by minor emission lines, which we neglected 
in our corrections of Section 3.2.2. We note that the magnitude 
of these offsets is of the order of the uncertainties in the SDSS 
AB calibration 5 . In the near-infrared, there is a tendency for the 
model to slightly overestimate the observed A^-band luminosity 
(by an amount corresponding to about 0.10 mag). This offset is 
well within the uncertainties of current stellar population synthesis 
models for this spectral region. We have checked that these offsets 
have a negligible influence on our results. 

To illustrate the quality of these fits, we show in Fig. 4 
three examples of the best-fit spectral energy distributions of 
galaxies spanning wide ranges in star formation and dust prop- 
erties in our sample, from quiescent, moderately dusty (IRAS 
F15 105+5959, top panel) to actively star-forming, highly dusty 
(IRAS F09253+1724, bottom panel). The middle panel shows the 
spectral energy distribution of a galaxy corresponding roughly to 



see http://www.sdss.Org/DR7/algorithms/fluxcal.html#sdss2ab 
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Figure 4. Best-fit models (in black) to the observed spectral energy distribu- 
tions (in red) of 3 galaxies spanning wide ranges in star formation and dust 
properties in the matched GAL£X-SDSS-2MASS-/RAS sample described 
in Section 2: from quiescent, moderately dusty (top panel) to actively star- 
forming, highly dusty (bottom panel). In each panel, the blue line shows the 
unattenuated stellar spectrum. The red squares show the observed broad- 
band luminosities (GALEX FUV, NUV; SDSS ugriz; 2MASS JHK a ; 
IRAS 12, 25, 60 and 100 /im) with their errors as vertical bars and the widths 



of the filters as horizontal bars. The fit residuals (L 
shown at the bottom of each panel. 



obs 



T mod 

A 



the median star formation and dust properties of the sample (IRAS 
F 15028+0820). 



4 THE RELATION BETWEEN STAR FORMATION 
ACTIVITY AND DUST CONTENT IN GALAXIES 

The method described in Section 3 above allows us to derive sta- 
tistical constraints on the star formation rate averaged over the last 
10 8 yr, if), the stellar mass, M», the total dust luminosity, L^ ot , 
the fraction of this contributed by dust in the diffuse interstellar 
medium, / M , and the total dust mass, Md, for the 3258 galaxies 
in the matched GALEX-SDSS-2MASS-IRAS sample described in 
Section 2. 

In Fig. 5, we show the resulting relation between dust mass 
M<j and star formation rate tp. The grey contour shows the distri- 
bution for the full sample, while grey points show the distribution 
of a sub-sample of 1658 galaxies with highest-S/N photometry at 
ultraviolet, optical and infrared wavelengths. This high-S/N sub- 
sample includes only galaxies with relative photometric errors less 
than 2a larger than the sample mean in each band. Also shown in 
the lower-right corner of Fig. 5 are the median errors in tp and Md, 
derived from the likelihood distributions of these parameters for all 
the galaxies in the sample. 




-3-2-10 1 2 3 
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Figure 5. Median-likelihood estimate of the dust mass, Md, versus that of 
the star formation rate averaged over the last 10 s yr, tp, for the matched 
GALEX-SDSS-2MASS-IRAS sample described in Section 2 (in grey) and 
the SINGS galaxies (in orange). The grey contour shows the distribution of 
the full matched sample of 3258 galaxies, while grey points show the distri- 
bution of the sub-sample of 1658 galaxies with highest-S/N photometry in 
all bands (see text). The median error bars for each sample are indicated in 
the lower-right corner. These correspond to the median 16th-84th percentile 
range of the likelihood distributions in Md and tp. The black line shows a 
linear fit to the grey points, computed as the bisector of the least-square 
regression lines on each axis. 



The correlation between dust mass and star formation rate is 
remarkably tight in Fig. 5, spanning 4 orders of magnitude in both 
tp and Md. We perform a linear fit to the grey points (i.e., the 
highest-S/N subsample) by computing the bisector of two least- 
square regression lines (Isobe et al. 1990): one in Md as a function 
of tp, and one in tp as a function of Md. This gives (black line in 
Fig. 5) 



M d = (1.28 ± 0.02) x 10 7 0>/Mq yr' 



-lll.lliu.01 



M, 







(9) 



The uncertainties quoted in this expression are formal ones derived 
by taking into account the confidence range in each tp and Md mea- 
surement. These uncertainties are much smaller than the intrinsic 
scatter of the relation (about 0.5 dex in Md). Equation (9) can thus 
be used to estimate the 'typical' dust mass in a galaxy, based on the 
star formation rate. To our knowledge it is the first time that such 
an expression is calibrated for a large sample of galaxies. 

We have investigated the extent to which the SDSS, GALEX 
and IRAS selections of our sample may introduce a bias in the rela- 
tion between star formation rate and dust mass derived from Fig. 5. 
We used for this the library of stochastic models described in Sec- 
tion 3.2.1. Since these models are normalised to total stellar mass 
(da Cunha et al. 2008), we assigned a random stellar mass (and 
scaled dust mass and star formation rate) to each model in the li- 
brary. We drew the stellar masses uniformly in log(M«/MQ) be- 
tween 8.5 and 1 1.5, to be consistent with the distribution of galaxy 
stellar masses in our sample. For each model in the library, we com- 
puted the expected GALEX, SDSS and IRAS magnitudes in differ- 
ent redshift bins from z — to z — 0.20. Then, we applied the 
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same selection criteria as used for our observed sample to the model 
library. These selection criteria introduce minimum detectable stel- 
lar mass, dust mass and star formation rate. We find that, for exam- 
ple, at the low-redshift end of our sample (z = 0.0025), the IRAS 
flux limit tends to exclude galaxies with log(Md/ Mq) < 4.5. 
At the typical redshift of our sample, z = 0.05, the combina- 
tion of SDSS and IRAS selections sets a minimum detectable dust 
mass of about 10 7 Mq and a minimum star formation rate of 
about 0.01 Mq yr _1 (it also introduces serious incompleteness 
for galaxies with star formation rates less than 1 Mq yr~ ). We 
have checked that these selection effects do not affect significantly 
the empirical relation between the star formation rate and dust mass 
derived from Fig. 5, which is dominated by more massive galaxies 
(including only galaxies with ip > 1 Mq yr _1 and Md > 10 7 Mq 
leaves the derived slope unchanged in equation 9). 

In Fig. 5, we also plot for comparison the dust masses and star 
formation rates derived by da Cunha et al. (2008) for the SINGS 
galaxies (orange symbols). These nearby galaxies follow the same 
relation as the galaxies in the matched GAL£X-SDSS-2MASS- 
IRAS sample studied here and extend to slightly lower dust masses 
and star formation rates 6 . We note that, at fixed star formation rate, 
the SINGS galaxies tend to have slightly larger dust masses than 
the galaxies in our sample. This is likely to result from the differ- 
ent selection criteria of the SINGS sample. The typical error bars 
in Md and ip are smaller for the SINGS galaxies than for the sam- 
ple studied here, because a wider collection of observational con- 
straints (especially in the mid-infrared) were available to da Cunha 
et al. (2008). 

In an attempt to understand which observations set the main 
constraints on dust mass, we have investigated how the median- 
likelihood estimates of Md correlate with a wide range of galaxy 
colours. We find that M d correlates most strongly with the 
F u wo /Fj> colour, where F„ 100 is the flux density in the IRAS 100- 
(im band and FJ that in the SDSS g band. The Spearman rank 
coefficient for this correlation is rs = 0.56, indicating a posi- 
tive correlation at the 23a level for this sample size (we note that 
the correlation is not improved by the inclusion of an ultraviolet 
band such as GALEX FUV). The reason for this correlation is that 
Fj 00 /FJ 1 traces the primary contributor to Md: cold dust in the 
diffuse ISM (Mq SM in eq. 7). This component is also mainly re- 
sponsible for the attenuation of the emission from stars older than 
10 7 yr, which dominate the SDSS g-band light, and it is the main 
contributor to the emission at 100 /im (a more detailed discussion 
of the set of observables necessary to constrain M d can be found in 
Section 3.2.2 of da Cunha et al. 2008). 

In Fig. 6, we explore the relations between specific star forma- 
tion rate tps and three particularly interesting physical properties of 
the galaxies in our sample: the dust-to-stellar mass ratio M d /M», 
the ratio of dust mass to star formation rate Md/ip, and the frac- 
tion / M of the total infrared luminosity L d ot contributed by dust in 
the diffuse ISM. As in Fig. 5, grey contours show the relations for 
the full sample of 3258 galaxies, while individual points show the 
relations for the 1658 galaxies with highest-S/N photometry. The 
top panels of Fig. 6 indicate that M d /M,, M d /ip and / M are all 
strongly correlated with : the Spearman rank coefficients for the 
full sample are rs = 0.84, —0.60 and —0.66, respectively, indi- 



° da Cunha et al. (2008) have checked that the dust masses of the SINGS 
galaxies derived using the model described in Section 3 above are typically 
within 50 per cent of those estimated by Draine et al. (2007) using a more 
sophisticated physical dust model. 



eating that the correlations are significant at more than 20ct level 
for this sample size (similar results are obtained when using only 
the high-S/N subsample). It is important to check that stellar mass 
is not the main driver for these strong correlations. To verify this, 
in the bottom panels of Fig. 6, we show differences in the same 
quantities as in the top panels for pairs of galaxies closely matched 
in stellar mass (we have included all possible galaxy pairs for any 
stellar mass M* in the full sample). Specifically, for each galaxy 
pair, we plot the difference in specific star formation rate between 
the two galaxies, [A log^s /yr -1 )], against the difference in dust- 
mass to stellar-mass ratio, [A log(Md /M* )], the difference in ratio 
of dust mass to star formation rate, [A log(M d /i/')], and the differ- 
ence in fraction of total infrared luminosity contributed by the dif- 
fuse ISM [A/ M ]. The fact that the strong correlations subsist from 
the top to the bottom panels of Fig. 6 demonstrates that stellar mass 
is not the main driver of these correlations. 

It is of interest to check the extent to which the properties of 
galaxies in Fig. 6 actually depend on stellar mass. Several stud- 
ies have shown that the star formation activity of a galaxy tends 
to decrease with increasing stellar mass (e.g., Brinchmann et al. 
2004). In Fig. 7, we plot M d /M*, M d /ip and / M as a function of 
ips for 3 different stellar-mass ranges chosen to contain roughly 
similar numbers of galaxies. The distributions in the various quan- 
tities on the j/-axis and their median values are displayed in the 
right-hand panels. Fig. 7 shows that the relation between Md/M* 
and ips depends little on stellar mass, with only a slight tendency 
for the less massive galaxies to have somewhat higher Mt/M*. 
Also, the quantity Md/ip, which traces the dust-to-gas ratio, tends 
to increase slightly with stellar mass. The quantity displaying the 
strongest dependence on stellar mass in Fig. 7 is the fraction / M of 
the total infrared luminosity contributed by dust in the diffuse ISM: 
the median value of this parameter increases from 0.5 for low-mass 
galaxies to 0.7 for high-mass ones. We note that the 3 stellar-mass 
bins considered here span only a small dynamic range in M* , which 
may not be enough to reveal the full dependence of the considered 
properties on stellar mass. It is tempting to interpret the relations 
of Figs. 6 and 7 as evolutionary sequences, where the build-up of 
stellar mass would be accompanied by gas consumption and dust 
enrichment. We return to this point Section 5.4 below. 

We have also checked that the non-detection in the 12- and 
25-/im IRAS bands of most galaxies in the sample (Section 2.4) 
has a negligible influence on the results of Fig. 6. This is illustrated 
by the blue triangles in the top panels of the figure, which mark 
galaxies with complete IRAS information in the 12-, 25-, 60- and 
100-/im bands. As expected, the median errors in the derived pa- 
rameters of these galaxies are typically smaller than those for the 
full sample (as indicated in the upper-right corner of each diagram). 
Fig. 6 shows that galaxies with complete IRAS information follow 
the same trends as the rest of the sample. Thus, the availability of 
only the 60- and 100-/im IRAS flux densities for most galaxies in 
the sample does not appear to systematically bias the results. 

The strength of the correlations between tps and Ma/M», 
Md/ip and / M in Fig. 6 suggests that the specific star formation 
rate is a fundamental diagnostic of the ISM properties in galaxies. 
The left-hand panels of Fig. 6 show that, for example, the galaxies 
with highest specific star formation rate are also the most dust-rich. 
The reason for this can be best understood with the help of the mid- 
dle panels. In these diagrams, Md/ip may be regarded as a proxy 
for the dust-to-gas ratio, since ip is tightly related to the gas mass 
(H I + H 2 ) by virtue of the Schmidt-Kennicutt law (Schmidt 1959; 
Kennicutt 1998a). The inverse correlation between M d /ip and tps 
suggests that the dust-to-gas ratio is larger in galaxies with low spe- 



Star formation activity and dust content in galaxies 9 




-4 -2 2 4 -4 -2 2 4 -4 -2 2 4 

A log ( fi I y~' ) A log ( f s / yr H ) A log ( fs / If** ) 



Figure 6. 7bp: Median-likelihood estimates of 3 galaxy properties against specific star formation rate, -0s . Left panel: ratio of dust mass to stellar mass, 
Af d / M„. Middle panel: ratio of dust mass to star formation rate, M d /tp (which may be used as a proxy for the dust-to-gas ratio; see Section 4). Right panel: 
fraction of total infrared luminosity contributed by dust in the ambient ISM, / M . In each panel, the grey contour shows the distribution of the full matched 
GALEX-SDSS-2MASS-/RAS sample described in Section 2, while the points show the distribution of the sub-sample of 1658 galaxies with highest-S/N 
photometry in all bands. The blue triangles show galaxies with complete IRAS information in the 12-, 25-, 60- and 100-^tm bands. The error bars represent 
the median confidence ranges in each parameter. Bottom: differences in the same quantities as in the top panels for pairs of galaxies closely matched in stellar 
mass (see text for detail). The Spearman rank coefficient rg is indicated in each panel. 



cific star formation rate, which have presumably exhausted their 
gas reservoir. In contrast, galaxies with high specific star formation 
rate and hence large gas reservoirs tend to have lower dust-to-gas 
ratio [we note that supernovae and young asymptotic-giant-branch 
(AGB) stars will produce large amounts of new dust in these galax- 
ies; see e.g. Dwek 1998]. The left-hand panels of Fig. 6 show that, 
as a result, a young, actively star-forming galaxy with low dust-to- 
gas ratio may still be highly dusty (in the sense of a high Md / M* ) 
because it contains large amounts of interstellar gas (see also Sec- 
tion 5.4 below). Finally, the right-hand panels of Fig. 6 show that 
stellar birth clouds are the main contributors to dust heating in ac- 
tively star-forming galaxies, while in more quiescent galaxies, the 
bulk of dust emission arises from the heating of dust by older stars 
in the diffuse interstellar medium. This trend was previously noted 
for SINGS galaxies by da Cunha et al. (2008). 

We may attempt to combine the above estimates of the dust 
mass Md for the galaxies in our sample with gas-mass estimates 
derived from the star formation rate ip via the Schmidt-Kennicutt 
law, to constrain the dust-to-gas ratio. The Schmidt-Kennicutt takes 
the form (Kennicutt 1998b) 

Espr = (2.5 ± 0.7) x 10~ 4 E H ' 4±0 - 15 , (10) 
where Esfr is the surface density of star formation, expressed in 



Mq yr _1 kpc -2 , and Eh is the surface mass density of H1+H2 
gas, expressed in Mq pc -2 . This formula was derived using 
a Salpeter (1955) IMF with lower and upper cut-offs 0.1 and 
100 Mq. We adjust the scaling coefficient by a factor of 1.6 down 
to account for the fact that the da Cunha et al. (2008) model used 
here relies on a Chabrier (2003) IMF with same cut-offs. We com- 
pute Esfr for the galaxies in the high-S/N subsample studied 
above by dividing our estimates of ip by the area defined by the 
minor and major axes of the r-band isophote at 25 mag arcsec -2 
(see Section 5.2). We then compute Eh and hence Mh using equa- 
tion (10). 

In Fig. 8, we plot the dust-to-gas ratio M^/Mh derived in this 
way against the stellar mass M* and the dust attenuation optical 
depth in the diffuse ISM firry for the galaxies in the high-S/N sub- 
sample. The median errors in M^/Mu are quite large. Yet, this ra- 
tio appears to correlate with both M* (rs = 0.33; 13er significance 
level) and fifv (rs = 0.50; 20a significance level). The correla- 
tion with M, is consistent with the tight correlation between stel- 
lar mass and gas-phase metallicity for SDSS star-forming galax- 
ies (e.g. Tremonti et al. 2004). In fact, the strong correlation be- 
tween dust-to-gas ratio and gas-phase metallicity has been noted 
in several previous studies (e.g. Issa et al. 1990; Schmidt & Boiler 
1993; Lisenfeld & Ferrara 1998). We note that the relation between 
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Figure 7. Same relations as in the top panels of Fig. 6 plotted in 3 stellar-mass ranges: log(M*/MQ) < 10.2 (blue); 10.2 ^ log(M» / Mq) < 10.6 
(orange); log(M* / Mq) > 10.6 (red). The right-hand panels show the corresponding distributions of the dust-to-stellar mass ratio Mj/M*, the ratio of dust 
mass to star formation rate M^/tp and the fraction / M of total infrared luminosity contributed by dust in the ambient ISM, along with their median values. The 
contours are the same as in the top panels of Fig. 6. 



Md /Mh and fifv in Fig. 8 has too much scatter for the dust opti- 
cal depth in the diffuse ISM to be used as a reliable proxy for the 
dust-to-gas ratio in galaxies. 



5 DISCUSSION 

In this section, we briefly mention model uncertainties and discuss 
potential biases introduced by inclination effects and the presence 
of obscured AGNs in the results derived in the previous sections. 
We also compare these results with the predictions of models of 
chemical and dust evolution to illustrate the implications of our 
study for the evolution of star-forming galaxies. 

5.1 Model uncertainties 

The da Cunha et al. (2008) model used in Section 3 to derive 
the physical properties of galaxies in the matched GALEX-SDSS- 
2MASS-IRAS sample relies on a combination of the latest version 
of the Bruzual & Chariot (2003) population synthesis code and the 



simple two-component dust model of Chariot & Fall (2000). Re- 
cently, Conroy et al. (2009) have investigated the uncertainties in- 
herent to such models, in particular, those arising from poorly un- 
derstood phases of stellar evolution (blue horizontal branch stars, 
blue straggler stars, thermally pulsing AGB stars), the IMF and the 
properties of dust in the interstellar medium. 7 Uncertainties linked 
to the IMF will affect our analysis in the same way as pointed out 
by these authors. In our case, the uncertainties linked to stellar evo- 
lution are somewhat minimized by the fact that blue-straggler and 



7 We note that Conroy et al. (2009) inappropriately refer to the Chariot & 
Fall (2000) model as a 'uniform screen of dust'. In reality, by design, the 
effective attenuation curve in the Chariot & Fall (2000) model reflects the 
probability density for the absorption of photons emitted in all directions 
by stars in all locations within a galaxy (see their eq. [18]). The shape of 
this curve - which is the quantity primarily constrained by observations - is 
controlled by the combination of the optical properties and the spatial dis- 
tribution of the dust. Allowing for variations in this shape therefore allows 
one to explore variations in the spatial distribution of the dust in a galaxy 
(as exemplified in section 4 of Chariot & Fall 2000). 
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Figure 8. Top: ratio of dust mass to gas mass, XI \\ = Mjj i + Mn 2 , 
plotted against the stellar mass M* for the 1658 galaxies with highest-S/N 
photometry in the matched GALEX-SDSS-2MASS-IRAS sample described 
in Section 2. Bottom: AI^/Mn ratio plotted against the dust attenuation 
optical depth in the diffuse ISM, firy, for the same sample. The median 
error bars are indicated in the lower-right corner of each panel. 



horizontal-branch stars have a negligible contribution to the emis- 
sion from star-forming galaxies. Also the version of the Bruzual & 
Chariot (2003) models used here rely on the improved, empirically 
calibrated models of Marigo & Girardi (2007) to better account 
for the contribution by thermally pulsing AGB stars to the near- 
infrared emission from galaxies (see Bruzual 2007). We refer the 
reader to the original studies of Chariot & Fall (2000) and da Cunha 
et al. (2008) for discussions of the uncertainties in the dust model 
(e.g., spatial distribution and optical properties of dust grains; typ- 
ical lifetime of giant molecular clouds) and their influence on the 
ultraviolet and infrared emission from galaxies. 



5.2 The effect of inclination 

The model used in Section 3 to derive the physical properties 
of galaxies in the matched GAL£X-SDSS-2MASS-//M5 sample 
is limited to angle-averaged spectral properties. However, the ob- 
served fluxes and colours of galaxies are observed to depend on 
inclination, especially for spiral galaxies (e.g., Mailer et al. 2009). 
SDSS galaxies have been morphologically classified using the con- 
centration index C, defined as the ratio of the radii enclosing 90 per 
cent and 50 per cent of the r-band luminosity (e.g. Shimasaku et al. 
2001; Strateva et al. 2001). We find that 64 per cent of the galaxies 
in our sample have concentration indices typical of late-type galax- 
ies (i.e., C < 2.6). Thus, we must check that inclination effects do 
not introduce any strong bias in the results of Section 4. 

We use the ratio b/a of minor to major axes of the SDSS r- 
band isophote at 25 mag arcsec -2 as a proxy for disc inclination. 
This is justified by the way in which galaxy ultraviolet and opti- 
cal colours correlate with b/a. For example, we find that colours 
such as g — r and NUV — r correlate significantly with b/a, with 
the smaller b/a (presumably tracing higher disc inclinations) corre- 
sponding to the redder ultraviolet/optical colours. In Fig. 9, we plot 
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Figure 9. Galaxy properties plotted against ratio of minor to major axes 
b/a (used as a proxy for inclination) for the 1658 galaxies with highest-S/N 
photometry in the matched GAL£X-SDSS-2MASS-/R45 sample described 
in Section 2: (a) apparent r-band magnitude, r; (b) stellar mass, M*; (c) 
star formation rate averaged over the last 10 s yr, i/>; (d) total infrared lumi- 
nosity, L^ ot ; (e) dust mass, Md ; (f) fraction of L^ ot contributed by dust in 
the diffuse ISM, fa ■ Red squares and error bars indicate the median values 
and 16th-84th percentile ranges in bins of b/a. For reference, the horizon- 
tal solid and dashed grey lines indicate the median value and 16th-84th 
percentile range for the whole sample (i.e. including all inclinations). The 
median observational error in r is indicated in the lower-left corner of panel 
(a). In each panel from (b) to (f), the median confidence interval in the like- 
lihood estimates of the considered parameter is indicated in the lower-left 
corner. 



several properties of the galaxies in the high-S/N subsample (Sec- 
tion 4) against the axis ratio b/a: apparent r-band magnitude, r; 
stellar mass, M* ; star formation rate averaged over the last 10 s yr, 
ip; total infrared luminosity 8 , L^ ot ; dust mass, Ma; and fraction of 
Ll ot contributed by dust in the diffuse ISM, / M . 

We first highlight a selection effect: at low b/a (i.e. high 
inclination), the sample includes almost no galaxy fainter than 
r ~ 16 mag (Fig. 9a). This is because of the combined require- 
ments of infrared detection by IRAS and ultraviolet detection by 
GALEX: edge-on galaxies in our sample must be optically bright 



8 We note that the values of L^ ot derived from our multi- wavelength anal- 
ysis are typically 0.2 dex larger than those computed using the simple em- 
pirical approximation of eqs. (l)-(2). 
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enough that the limited relative dust content required for some ul- 
traviolet photons to escape in the plane of the disc be also sufficient 
to warrant detection in the infrared. It is important to note that our 
restriction to the high-S/N subsample in Fig. 9a explains the lack 
of galaxies at the faintest magnitudes near the r = 17.77 selection 
cut-off even at low inclinations. 

Fig. 9b shows that the median-likelihood estimates of stellar 
mass do not vary significantly with axis ratio, implying that incli- 
nation has a negligible influence on M* constraints. This is not sur- 
prising, since constraints on M* are dominated by the near-infrared 
emission, which is little sensitive to dust attenuation. Estimates of 
the star formation rate exhibit a stronger dependence on inclination 
(Fig. 9c): the typical ip of face-on galaxies (b/a « 1) is a factor 
of about 3 larger than that of edge-on galaxies. This difference is 
significant when compared to the median confidence interval in ip 
estimates (indicated on the figure). It is consistent with the selec- 
tion bias identified in Fig. 9a against very dusty edge-on galaxies. 
This effect is also responsible for the lower dust masses and dust 
luminosities derived for edge-on compared to face-one galaxies in 
Figs. 9d and 9e. Finally, Fig. 9f shows the contribution / M to the to- 
tal infrared luminosity by dust in the diffuse ISM is systematically 
higher in edge-on galaxies than in face-on ones. This is likely to 
result from the lower typical specific star formation rate of edge-on 
galaxies (Figs. 9b and 9c; see the correlation between / M and in 
the right-most panels of Fig. 6). 

We emphasise that the effects identified above are weak com- 
pared to the typical dispersion in the estimated parameters. Thus, 
inclination effects can account only for a minor fraction of the in- 
trinsic scatter in the relations studies in Section 4. 

5.3 Contamination by AGN hosts 

In Section 3, we have interpreted the spectral energy distributions 
of galaxies in the matched GALEX-SDSS-2MASS-IRAS sample 
using the model of da Cunha et al. (2008). An important limita- 
tion of this model is that it does not include the potential contri- 
bution by an AGN to the infrared emission of a galaxy (e.g. de 
Grijp et al. 1985). For this reason, in Section 2, we have excluded 
from our sample all potential AGN hosts on the basis of their 
optical-line emission. However, optically-thick AGNs would not 
be eliminated using this method. This question is a fair concern, 
as the fraction of AGN hosts is known to increase with the total 
infrared luminosity (Veilleux et al. 1995; Cao et al. 2006). In addi- 
tion, evidence for dust-enshrouded AGN activity has been found in 
some ULIRGs and LIRGs of the IRAS catalogue (Dudley & Wynn- 
Williams 1997; Imanishi & Dudley 2000). In an infrared-selected 
sample of high-redshift (z ~ 0.8) galaxies with infrared luminosi- 
ties similar to those of the galaxies in our sample, Elbaz et al. 
(2002) find 12 ± 5 per cent (5/41) of AGN hosts based on avail- 
able observations with Chandra X-ray Observatory (Brandt et al. 
2001; Hornschemeier et al. 2001). It is conceivable, therefore, that 
a similar fraction of galaxies in our sample host an obscured AGN. 

Hot dust in the torus surrounding an obscured AGN is ex- 
pected to radiate mostly in the mid-infrared. This should have only 
a minor influence on our analysis, since most galaxies in our sam- 
ple are not detected in the IRAS 12- and 25-/im bands (Fig. 6). 
However, some studies suggest that AGNs may be also contribute 
significantly to the far-infrared emission from galaxies. For exam- 
ple, in the analysis of a matched sample of SDSS-DR2 and IRAS 
galaxies, Pasquali et al. (2005) find that known AGN hosts exhibit 
a significant excess of 60- and 100-/im emission, amounting to a 
typical excess infrared luminosity of 0.18 dex. Yet, it remains un- 



clear whether this excess emission arises from (cold) dust heated 
by the AGN or by a population of young stars not detected in the 
optical SDSS spectra. In another recent study, Salim et al. (2009) 
analyse the 24 fim emission of a sample of optically-selected galax- 
ies observed with the Spitzer Space Telescope (Werner et al. 2004) 
out to redshift z = 1.4. They conclude that AGNs identified in the 
X-rays and the optical do not contribute significantly to the 24 ^im 
emission in their sample. Moreover, optically selected AGNs do 
not present any evidence for an excess infrared luminosity, while 
obscured AGNs detected in the X-rays can contribute only up to 
50 per cent of the total infrared luminosity. 

We have checked that the possible excess infrared luminosity 
arising from an obscured AGN in these previous studies is within 
the typical confidence range of estimates for the galaxies in 
our sample. We conclude that the potential contamination of the 
60- and 100-/im IRAS flux densities by obscured AGNs would not 
alter significantly our conclusions. 

5.4 Comparison with models of chemical and dust evolution 

In this section, we compare the physical properties derived in Sec- 
tion 4 for the galaxies in the matched GALEX-SDSS-2MASS-IRAS 
sample with the predictions of chemical evolution models includ- 
ing the formation and destruction of dust grains in the ISM. We 
focus on the recent models of the chemical and dust evolution of 
galaxies in a wide range of star formation histories by Calura et al. 
(2008). These models, which are based on the formalism developed 
by Dwek (1998), include the production of dust in the cold en- 
velopes of low- and intermediate-mass stars during the AGB phase 
and in the expanding ejecta of type-II and type-la supernovae. Dust 
grains are destroyed mainly by shock waves caused by supernovae 
explosions. The Calura et al. (2008) models allow one to compute 
the evolution of the stellar and dust masses of galaxies as a function 
of the star formation history. 

We are particularly interested here in the Calura et al. (2008) 
models for the solar neighbourhood and for a dwarf-irregular 
galaxy with continuous star formation. The solar-neighbourhood 
model aims at reproducing the properties of the Galactic-disc com- 
ponent in a 2 kpc-wide ring located at 8 kpc from the Galactic cen- 
tre. In this model, which assumes gas infall but no outflow, the 
Galactic disc is formed in two main episodes. The exact expression 
of the star formation history, taken from Chiappini et al. (1997), is 
essentially a function of the gas surface mass density with a crit- 
ical threshold for star formation. In the model of dwarf-irregular 
galaxy, the star formation rate is expressed as a Schmidt-Kennicutt 
law and increases continuously with time. This model assumes gas 
infall and galactic winds. We refer the reader to the original paper 
of Calura et al. (2008) for more detail on these models. 

In Fig. 10, we compare the median-likelihood estimates of 
Md/M, and Ma/ip as a function of tps for the galaxies in the 
high-S/N subsample of Section 4 with the predictions of the solar- 
neighbourhood model (red filled squares) and the dwarf-irregular 
galaxy model (blue filled circles) of Calura et al. (2008), at ages be- 
tween 1 and 12 Gyr. At late ages, both models coincide remarkably 
well with the range of derived physical parameters of the GALEX- 
SDSS-2MASS-WMS galaxies. The solar-neighbourhood model of 
Calura et al. (2008) does not include any bulge component. To lo- 
cate the Milky Way in Fig. 10, we add to this model a bulge stellar 
mass of 1/3 that of the disc. We assume that this bulge compo- 
nent has negligible contributions to M d and ips (consistent with 
the chemical evolution model of a 12-Gyr old spheroidal galaxy by 
Pipino et al. 2005). The resulting Milky-Way model has a total stel- 
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Figure 10. Comparison of the physical properties derived in Section 4 for the high-S/N galaxies in the matched GALEX-SDSS-2MASS-IRAS sample described 
in Section 2 (grey points) with the theoretical predictions of the Calura et al. (2008) models of chemical and dust evolution. Left: dust-to-stellar mass ratio 
as a function of specific star formation rate i/)g. Right: ratio of dust mass to star formation rate M^/ip as a function of ^>g. The evolution of the 
solar-neighbourhood model of Calura et al. (2008) is shown in red, and that of the dwarf-irregular galaxy model in blue. Symbols indicate the locations of 
these models at ages 1, 3, 6 and 12Gyr. The orange arrow shows the effect of adding a bulge component to the 12-Gyr old solar-neighbourhood model to 
mimic the Milky- Way properties (see text for detail). 



lar mass M, — 4 x 10 10 Mq, a star formation rate tj) = 2.3 Mq 
yr -1 and a dust mass Md = 4.5 x 10 7 Mq (head of the orange 
arrow in Fig. 10). These parameters are typical of moderately star- 
forming and moderately dusty galaxies in our sample. 



A valuable feature of the models shown in Fig. 10 is that 
they provide a framework to interpret the time evolution of the 
ISM properties of galaxies. Both the solar-neighbourhood and the 
dwarf-irregular galaxy models predict that, overall, Md /M* should 
drop and M d /V> should rise as the specific star formation rate 
ips declines. These trends, and the good general agreement of the 
models with the derived physical properties of the galaxies in our 
sample in Fig. 10, suggest that the relations identified between 
M d /M*, Md/ip, / p and ips in Section 4 (Fig. 6) may be the result, 
at least in part, of an evolutionary sequence. In detail, the interpre- 
tation of this sequence depends on several key assumptions of the 
chemical evolution models, such as the star formation history, the 
influence of gas infall and outflow, and the dust production and de- 
struction rates. In the Calura et al. (2008) models of Fig. 10 (which 
both include gas infall, the dwarf irregular model also including 
gas outflow), a young galaxy has an ISM characterized by large 
amounts of dust-poor gas and forms stars at a very high rate. As 
the galaxy continues to consume gas into star formation, the stellar 
mass rises and the star formation rate declines, causing ips to drop. 
The dust mass increases steeply during the first billion years of evo- 
lution, when the star formation rate is highest and large amounts of 
dust are produced in supernova ejecta and in the envelopes of AGB 
stars. As star formation drops, the dust mass Md also drops, be- 
cause the production of new dust grains cannot balance any more 
their destruction in the ISM. This general picture, which is consis- 
tent with our results, is also in broad agreement with the predictions 
of simple closed-box chemical evolution models (e.g., Dwek et al. 
2000; Edmunds 2001). 



6 SUMMARY AND CONCLUSION 

In this paper, we have assembled a sample of 3258 local galaxies 
with photometric observations from at ultraviolet, optical and in- 
frared wavelengths. These galaxies were primarily selected from 
the SDSS DR6 spectroscopic sample cross-correlated with IRAS 
all-sky catalogues (which are flux-limited at 60 /im). The optical 
and infrared observations were supplemented with matched ultravi- 
olet and near-infrared data from GALEX and 2MASS, respectively. 
The size of this sample represents a significant improvement over 
previous multi- wavelength studies of galaxies. We have used the 
model of da Cunha et al. (2008) to interpret the observed spectral 
energy distributions of the galaxies in this sample in terms of sta- 
tistical constraints on the star formation rate, stellar mass, dust at- 
tenuation, dust luminosity, fraction of dust luminosity contributed 
by the diffuse ISM, and dust mass. 

We have focused on a subsample of 1658 galaxies with 
highest-S/N photometry to investigate several significant correla- 
tions between various derived physical properties of galaxies. In 
particular, we find that the star formation rate averaged over the last 
10 s yr, tp, correlates remarkably well with galaxy dust mass Md 
over 4 orders of magnitude in both quantities. The simple empiri- 
cal recipe M d = (1.28±0.02) x 10 7 0>/Mq yr" 1 ) 1 11±0 01 Mq 
may be used to roughly estimate the total dust mass of a galaxy as 
a function of the star formation rate. We also find that the dust-to- 
stellar mass ratio Md /M* , the ratio of dust mass to star formation 
rate Md / tp and the fraction / M of dust luminosity contributed 
by the diffuse ISM correlate strongly with the specific star forma- 
tion rate ips for the galaxies in the high-S/N subsample. Some of 
these trends had already been anticipated in the analysis of a much 
smaller sample of 66 SINGS galaxies by da Cunha et al. (2008). 

To investigate the origin of these correlations between vari- 
ous physical properties of galaxies, we have compared our results 
with the predictions of recent models of chemical and dust evolu- 
tion of galaxies with different star formation histories by Calura 
et al. (2008). We conclude from this comparison that the relations 
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between M^/M,, M^/tp, / M and tps could arise, at least in part, 
from an evolutionary sequence. As galaxies form stars, their ISM 
becomes enriched in dust, while the drop in gas supply makes the 
specific star formation rate decrease. Interestingly, as a result, a 
young, actively star-forming galaxy with low dust-to-gas ratio may 
still be highly dusty (in the sense of a high M d /M») because it 
contains large amounts of interstellar gas (Fig. 6). This may be im- 
portant for the interpretation of the infrared emission from young, 
gas-rich star-forming galaxies at high redshift. 

The results presented in this paper should be especially useful 
to improve the treatment of the ISM properties of galaxies in semi- 
analytic models of galaxy formation. More specifically, the results 
of Section 4 provide interesting constraints on the relation between 
star formation activity and dust content in such models. For exam- 
ple, the derived relation between Ma and ip (Fig. 5) can be used in 
combination of the Schmidt-Kennicutt law (Schmidt 1959; Kenni- 
cutt 1998a) to describe the evolution of the dust content of galaxies 
together with that of the stars and gas. Also, the correlation between 
tps and fn (Fig. 6) provides valuable clues on the way in which to 
implement the modelling of the different phases of the ISM in sim- 
ulated galaxies. 

Our study provides a local reference for comparison with 
the properties of galaxies at high redshifts. Observations with the 
Spitzer satellite have shown that infrared galaxies dominate the star 
formation activity of the universe at redshift z ~ 1 (Le Floc'h 
et al. 2005). Moreover, dust emission has been detected at sub- 
millimetre wavelengths out to z ~ 6, implying large dust masses 
and star formation rates of the order of several x 10 3 Mq yr _1 in 
(at least some) young galaxies (e.g., Chini & Kruegel 1994; Isaak 
et al. 2002; Blain et al. 2002; Bertoldi et al. 2003; Walter 2009). 
The production of such large amounts of dust on short timescales 
might require both efficient star formation and efficient grain con- 
densation in supernova remnants (e.g., Dunne et al. 2003; Morgan 
& Edmunds 2003; Maiolino et al. 2004). This could affect the re- 
lation between star formation rate and dust mass. Forthcoming ob- 
servations with the Herschel Space Telescope and the James Webb 
Space Telescope will enable the extension to high redshift of the 
type of analysis achieved in this paper. This can only lead to an im- 
proved understanding of the cosmic evolution of the dust content 
of galaxies. 



ACKNOWLEDGEMENTS 

We thank the referee, Charlie Conroy, for comments and sugges- 
tions which helped improve the quality of this paper. We are deeply 
grateful to David Schiminovich for matching the SDSS-/&4S sam- 
ple with the latest GALEX data release and for kindly providing us 
with these data. We thank Jarle Brinchmann for useful discussions 
and for providing us with the corrections to account for emission- 
line contamination in the optical SDSS bands. We also thank 
Loretta Dunne, David Elbaz, Vivienne Wild, Jakob Walcher and 
Brent Groves for useful discussions. We thank Francesco Calura 
and Antonio Pipino for providing us with the electronic data of their 
chemical evolution models. EdC was financed by the EU Marie 
Curie Research Training Network MAGPOP and by the EU Marie 
Curie ToK Grant MTKD-CT-2006-039965. CE was supported by 
the Marie Curie EARA-EST host fellowship while this work was 
carried out and acknowledges the IAP for hospitality. CE is partly 
supported by the Swiss Sunburst Fund. 

Funding for the SDSS has been provided by the Alfred P. 
Sloan Foundation, the Participating Institutions, the National Sci- 



ence Foundation, the US Department of Energy, the National Aero- 
nautics and Space Administration, the Japanese Monbukagakusho, 
the Max Planck Society, and the Higher Education Funding Coun- 
cil for England. The SDSS Web site is http://www.sdss.org. The 
SDSS is managed by the Astrophysical Research Consortium for 
the Participating Institutions. The Participating Institutions are the 
American Museum of Natural History, the Astrophysical Insti- 
tute Potsdam, the University of Basel, Cambridge University, Case 
Western Reserve University, the University of Chicago, Drexel 
University, Fermilab, the Institute for Advanced Study, the Japan 
Participation Group, Johns Hopkins University, the Joint Institute 
for Nuclear Astrophysics, the Kavli Institute for Particle Astro- 
physics and Cosmology, the Korean Scientist Group, the Chinese 
Academy of Sciences, Los Alamos National Laboratory, the Max 
Planck Institute for Astronomy, the Max Planck Institute for Astro- 
physics, New Mexico State University, Ohio State University, the 
University of Pittsburgh, the University of Portsmouth, Princeton 
University, the US Naval Observatory, and the University of Wash- 
ington. 



REFERENCES 

Adelman-McCarthy et al. 2008, ApJS, 175, 297 
Armus et al. 2007, ApJ, 656, 148 

Baldwin J. A., Phillips M. M., Terlevich R., 1981, PASP, 93, 5 
Beichman C. A., Neugebauer G, Habing H. J., Clegg P. E., 

Chester T. J., eds, 1988, Infrared astronomical satellite (IRAS) 

catalogs and atlases. Volume 1: Explanatory supplement Vol. 1 
Bertoldi F, Carilli C. L., Cox P., Fan X., Strauss M. A., Beelen 

A., Omont A., Zylka R., 2003, A&A, 406, L55 
Blain A. W., Smail I., Ivison R. J., Kneib J.-P, Frayer D. T, 2002, 

Physics Reports, 369, 111 
Blanton et al. 2003, AJ, 125, 2348 

Brandt W. N., Hornschemeier A. E., Alexander D. M., Garmire 
G. P., Schneider D. P., Broos P. S., Townsley L. K., Bautz M. W, 
Feigelson E. D., Griffiths R. E., 2001, AJ, 122, 1 
Brinchmann J., Chariot S., White S. D. M., Tremonti C, Kauff- 
mann G, Heckman T, Brinkmann J., 2004, MNRAS, 351, 1 151 
Bruzual G, 2007, (astro-ph/0703052) 
Bruzual G, Chariot S., 2003, MNRAS, 344, 1000 
Calura F, Pipino A., Matteucci F, 2008, A&A, 479, 669 
Cao C, Wu H., Wang J.-L., Hao C.-N., Deng Z.-G, Xia X.- 
Y, Zou Z.-L., 2006, Chinese Journal of Astronomy and Astro- 
physics, 6, 197 
Chabrier G, 2003, PASP, 115, 763 
Chariot S„ Fall S. M., 2000, ApJ, 539, 718 
Chiappini C, Matteucci F, Gratton R., 1997, ApJ, 477, 765 
Chini R., Kruegel E„ 1994, A&A, 288, L33 
Conroy C, White M., Gunn J. E., 2009, ArXiv e-prints 
da Cunha E., Chariot S., Elbaz D., 2008, MNRAS, 388, 1595 
de Grijp M. H. K., Miley G. K., Lub J„ de Jong T, 1985, Nature, 
314,240 

Dopita M. A., Groves B. A., Fischera J., Sutherland R. S., Tuffs 

R. J., Popescu C. C, Kewley L. J., Reuland M., Leitherer C, 

2005, ApJ, 619, 755 
Draine B. T, Dale D. A., Bendo G, et al. 2007, ApJ, 663, 866 
Dudley C. C, Wynn- Williams C. G, 1997, ApJ, 488, 720 
Dunne L., Eales S., Edmunds M., Ivison R., Alexander P., 

Clements D. L., 2000, MNRAS, 315, 115 
Dunne L., Eales S., Ivison R., Morgan H., Edmunds M., 2003, 

Nature, 424, 285 



Star formation activity and dust content in galaxies 15 



Dwek E., 1998, ApJ, 501, 643 

Dwek E., Fioc M., Varosi E, 2000, in D. Lemke, M. Stickel, & 
K. Wilke ed., ISO Survey of a Dusty Universe Vol. 548 of Lec- 
ture Notes in Physics, Berlin Springer Verlag, The Effect of Dust 
Evolution on the Spectral Energy Distribution of Galaxies, pp 
157-+ 

Edmunds M. G., 2001, MNRAS, 328, 223 

Elbaz D., Cesarsky C. J., Chanial P., Aussel H., Franceschini A., 

Fadda D., Chary R. R., 2002, A&A, 384, 848 
Helou G., Khan I. R., Malek L., Boehmer L., 1988, ApJS, 68, 151 
Hildebrand R. H., 1983, QJRAS, 24, 267 

Hopkins A. M., Connolly A. J., Haarsma D. B., Cram L. E., 2001, 
AJ, 122, 288 

Hornschemeier A. E., Brandt W. N., Garmire G. P., Schneider 
D. P., Barger A. J., Broos P. S., Cowie L. L., Townsley L. K., 
Bautz M. W., Burrows D. N., Chartas G., Feigelson E. D., Grif- 
fiths R. E., Lumb D., Nousek J. A., Ramsey L. W., Sargent 
W. L. W., 2001, ApJ, 554, 742 

Imanishi M., Dudley C. C, 2000, ApJ, 545, 701 

Isaak K. G., Priddey R. S., McMahon R. G., Omont A., Peroux 
C, Sharp R. G„ Withington S., 2002, MNRAS, 329, 149 

Isobe T., Feigelson E. D., Akritas M. G., Babu G. J., 1990, ApJ, 
364, 104 

Issa M. R., MacLaren I., Wolfendale A. W., 1990, A&A, 236, 237 
Kauffmann et al. 2003a, MNRAS, 341, 33 
Kauffmann et al. 2003b, MNRAS, 341, 54 
Kennicutt Jr. R. C, 1998a, ARA&A, 36, 189 
Kennicutt Jr. R. C, 1998b, ApJ, 498, 541 

Kennicutt Jr. R. C, Armus L., Bendo G., Calzetti D., Dale D. A., 
Draine B. T., Engelbracht C. W., Gordon K. D., et al. 2003, 
PASP, 115, 928 

Kewley L. J„ Jansen R. A., Geller M. J., 2005, PASP, 1 17, 227 
Kong X., Chariot S., Brinchmann J., Fall S. M., 2004, MNRAS, 
349, 769 

Le Floc'h E., Papovich C, Dole H., et al. 2005, ApJ, 632, 169 

Lisenfeld U., Ferrara A., 1998, ApJ, 496, 145 

Maiolino R., Schneider R., Oliva E., Bianchi S., Ferrara A., Man- 

nucci F, Pedani M., Roca Sogorb M., 2004, Nature, 431, 533 
Mailer A. H., Berlind A. A., Blanton M. R., Hogg D. W„ 2009, 

ApJ, 691, 394 
Marigo P., Girardi L„ 2007, A&A, 469, 239 
Martin et al. 2005, ApJ, 619, LI 

Morgan H. L., Edmunds M. G., 2003, MNRAS, 343, 427 
Morrissey et al. 2005, ApJ, 619, L7 

Moshir M., 1989, IRAS Faint Source Survey, Explanatory sup- 
plement version 1 and tape. Pasadena: Infrared Processing and 
Analysis Center, California Institute of Technology, 1989, edited 
by Moshir, M. 

Obric M., Ivezic Z., Best P. N, Lupton R. H, Tremonti C, Brinch- 
mann J., Agueros M. A., et al. 2006, MNRAS, 370, 1677 

Pasquali A., Kauffmann G., Heckman T. M., 2005, MNRAS, 361, 
1121 

Pipino A., Kawata D., Gibson B. K, Matteucci F, 2005, A&A, 
434, 553 

Rigopoulou et al. 1999, in Cox P., Kessler M., eds, The Universe 
as Seen by ISO Vol. 427 of ESA Special Publication, Ultralumi- 
nous IRAS galaxies as seen with ISO. pp 833-+ 

Salim S., Dickinson M., Rich R. M., Chariot S., Lee J. C, Schimi- 
novich D., Perez-Gonzalez P. G., Ashby M. L. N., Papovich C, 
Faber S. M., Ivison R. J., Frayer D. T, Walton J. M., Weiner 
B. J., Chary R.-R., Bundy K, Noeske K, Koekemoer A. M., 
2009, ArXiv e-prints 



SalpeterE. E., 1955, ApJ, 121, 161 
Saunders et al. 2000, MNRAS, 317, 55 

Schlegel D. J., Finkbeiner D. P., Davis M., 1998, ApJ, 500, 525 
Schmidt K.-H., Boiler T, 1993, Astronomische Nachrichten, 314, 
361 

Schmidt M., 1959, ApJ, 129, 243 
Seibert et al. 2005, ApJ, 619, L55 

Shimasaku K, Fukugita M., Doi M., Hamabe M., Ichikawa T, 
Okamura S., Sekiguchi M., Yasuda N., Brinkmann J., Csabai I., 
Ichikawa S., Ivezic Z., Kunszt P. Z., Schneider D. P., Szokoly 
G. P., Watanabe M., York D. G., 2001, AJ, 122, 1238 

Silva L., Granato G. L., Bressan A., Danese L., 1998, ApJ, 509, 
103 

Skrutskie et al. 2006, AJ, 131, 1163 

Soifer B. T, Sanders D. B., Madore B. F, Neugebauer G., Daniel- 
son G. E., Elias J. H., Lonsdale C. J., Rice W. L., 1987, ApJ, 320, 
238 

Strateva et al. 2001, AJ, 122, 1861 

Sullivan M., Mobasher B., Chan B., Cram L., Ellis R., Treyer M., 

Hopkins A., 2001, ApJ, 558, 72 
Tremonti et al. 2004, ApJ, 613, 898 

Veilleux S., Kim D.-C, Sanders D. B., 1999, ApJ, 522, 1 13 
Veilleux S., Kim D.-C, Sanders D. B., Mazzarella J. M., Soifer 

B.T, 1995, ApJS, 98, 171 
Walter F, 2009, ArXiv e-prints 
Wang B., Heckman T. M., 1996, ApJ, 457, 645 
Werner M. W., Roellig T. L., Low F. J., Rieke G. H., Rieke M., 

Hoffmann W. E, Young E., Houck J. R., 2004, ApJS, 154, 1 



